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In previous studies, we have explored the ansatz that the volume elements of the Bures metrics 
over quantum systems might serve as prior distributions, in analogy with the (classical) Bayesian 
role of the volume elements ("Jeffreys' priors") of Fisher information metrics. Continuing this work, 
we obtain exact Bures prior probabilities that the members of certain Zotti-dimensional subsets of 
the /j/feen-dimensional convex set of 4 x 4 density matrices are separable or classically correlated. 
The main analytical tools employed are symbohc integration and a formula of Dittmann (J. Phys. 
A 32, 2663 [1999]) for Bures metric tensors. This study complements an earlier one (J. Phys. A 32, 
5261 [1999]) in which numerical (randomization) — but not integration — methods were used to 
estimate Bures separability probabilities for unrestricted 4x4 and 6x6 density matrices. The exact 
values adduced here for pairs of quantum bits (qubits), typically, well exceed the estimate (~ .1) 
there, but this disparity may be attributable to our focus on special low-dimensional subsets. Quite 
remarkably, for the q ~ 1 and q ~ \ states inferred using the principle of maximum nonadditive 
(Tsallis) entropy, the Bures probabilities of separability are both equal to %/2 — 1. For the Werner 
qubit-qutnt and qutrit-qutnt states, the probabilities are vanishingly small, while in the qubit-qubit 
case it is j. 
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I. INTRODUCTION 



A. Background 



In a previous study Q , we exploited certain numerical methods to estimate the a priori probabihty — based on the 
volume element of the Bures metric — that, a member of the fifteen-dimensional convex set (i?) of 4 x 4 density 
matrices is separable (classically correlated), that is, expressible as a convex combination of tensor products of pairs 
of 2 X 2 density matrices. (Ensembles of separable states, as well as of bound entangled states can not be "distilled" 
to obtain pairs in singlet states for quantum information processing 0,^.) This Bures probability estimate ~ .1 was 
rather unstable in character Table 1], due, in part it appeared, to difficult-to-avoid "over-parameterizations" of R, 



as well as to the unavailability, in that context, of numerical integration methods. But now in sees. [I A, 11 B and II C 
below, we are able to report exact probabilities of separability by restricting consideration to certain low-dimcnsional 
subsets of R, for which symbolic integration can be performed. Then, in the subsequent body of the paper, we 
investigate analogous questions when R is replaced by the convex sets 9x9 and 6x6 density matrices. (We have 
also studied, using numerical methods, the Bures probability of separability of the two-party Gaussian states |^ (cf. 
0)-) 




Preliminarily though, we investigate in sec. IB certain relevant motivating issues, first having arisen in the context 
of the 3x3 density matrices. These quantum-theoretic entities belong to an eight-dimensional convex set (Q), which 
we parameterize in the manner, 

-I / V + z u ~ iw X — iy\ 

(1) 

The feasible range of the eight parameters — defined by the boundary of Q — is determined by the requirements 
imposed on density matrices, in general, that they be Hermitian, nonnegative definite (all eigenvalues nonnegative) , 
and have unit trace . 

Dittmann [|[ eq. (3.8)] (cf. j^]) has presented an "explicit" formula (one not requiring the computation of eigenvalues 
and eigenvectors) for the Bures metric ( (2[j^) over the 3x3 density matrices. It takes the form 

dBuresip, P + dpf = ^Trjdpdp + ^-1-3 (dp - pdp){dp - pdp) + -^-^M_{dp - p-'dp){dp - p-'dp)}. (2) 

If we implement this formula, using pg for p, we obtain an 8 x 8 matrix — the Bures metric tensor, which we will 
denote by g. 

It has been proposed ||l|,[l2|-|l6|] that the square root of the determinant of g, that is, Iffl^^^, which gives the volume 
element of the metric, be taken as a prior distribution (to speak in terms of the specific instance presently before us) 
over the 3x3 density matrices (cf. This ansatz is based on an analogy with Bayesian theory |18 1^, in which 



the volume element of the Fisher information [pO| matrix is used as a reparameterization-invariant prior, termed 
"Jeffreys' prior". 

Unfortunately, the ("brute force") computation of the determinant of such 8x8 symbolic matrices appears to 



exceed present capabilities |21 2^. In light of this limitation, we pursued a strategy of fixing (in particular, setting to 
zero) a certain number (four) of the eight parameters, thus, leading to an achieveable calculation. A similar course 
was followed in a brief exercise in eqs. (31), (32)], but using a quite different parmeterization of Q — one based on 
the expected values with respect to a set of four mutually unbiased (orthonormal) bases of three-dimensional Hilbert 
space JglQ. 

In ]l6|] we have reported exact results for the "Hall normalization constants" for the Bures volumes of the n-state 
quantum systems, n = 2, . . . , 6. These analyses utilized certain parameterizations (of Schur form) of the nx n density 
matrices [ p5[ , in which the eigenvalues and eigenvectors of these density matrices are explicitly given. It was established 
there |]l^, sec. II. B], among other things, that the Bures volume element for the 3x3 density matrices is, in fact, 
normalizable over Q, forming a probability distribution. Since it appears to be highly problematical to find an explicit 
transformation from this eigen-parameterization of Boya et al p5[ ] to that used in (|l]) , we can not conveniently utilize 
the results of | p6| for our specific purposes here. For the n-level systems, n > 3, the analogous task would appear 
to be even more challenging, since the appropriate parameterizations of SU{n) and the associated invariant (Haar) 



measures seem not yet to have been developed (cf. [P6|-p9[). 
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B. Two forms of conditional Bures priors for a /oiir-parameter i/iree-level system 



Since we aim to reexamine the specific findings in [^2[, we will thus cast our analyses specifically in terms of the 
parameterized form (|^). We have computed, using the formula (H), the 8x8 Bures metric tensor (g) associated 
with (|). Then — only subsequent to this computation — we set the four parameters s, t, u and v all equal to zero 
in 5, obtaining what we denote by g. (Actually, this "conditioning" on P can be performed immediately after the 
determination of the differential element dp, thus simplifying the further calculations in (H).) Then, 

1 

(3) 



64w(l — w) 2 iy'^ — 



z2)5(a;2 +y'^ + z'^ - {v - 2)2) ' 



which could be considered to constitute the (unnormalized) conditional Bures prior over P. 

On the other hand, if we ab initio nullify the same four parameters {s,t,u,v) in pg, we get the family of density 
matrices, defined over a four-dimensional convex subset (P) of Q, 



V + z 







12/ 



X — ly 
2-2v 
V- 



(4) 



which was the specific object of study in |l2| . 

Now, let us describe two ways in which an alternative to the presumptive conditional prior (^, that is. 



16f (1 — 7j) 2 (ti2 — — y2 — z^) : 



(5) 



has been derived. We obtain the outcome (||) if we either: (a) employ pp directly in Dittmann's formula (||), and 
generate the corresponding 4x4 metric tensor and compute the square root of its determinant (the procedure followed 
in [l^); or (b) extract from g the 4x4 submatrix with rows and columns associated with (the four now- nullified 
parameters) d, x, y and z, that is 



4(i;2 



z2) 



-y 



xy 



Ml 



V 



7 



(6) 



and calculate the square root of its determinant. The matrix (|^) is not exactly the same as (14) there — a result 
which was presumably also obtained by the use of pp in (|^. The four diagonal entries there are the negatives of the 
ones in (0), that is, we had there fl^, eq. (14)] 



4(t;2 



— y2 — 2:2) 



V 



y +z 



-y 

xy 



X +z 



\ 



J 



(7) 



In any case, the determinants of these two nonidentical 4x4 matrices are the same, so the substantive conclusions of 
[ p2| regarding Bures priors are unchanged. 

M. J. W. Hall has pointed out that the result (||) is, in fact, an eight-dimensional volume element rather that the 
four-dimensional one desired here. In addition, a referee has remarked there can be "two different sets of basis one- 
forms that are used to compute the volume element. This happens, for example, in SU{n) when one uses A~^dA as 
the matrix of left invariant one-forms. There exist invariant forms in this matrix. One must choose an independent 
set. This set is thus not unique" . 

It is interesting to compare the form of (^ with the Bures metric tensor for the 2x2 systems eq. (4)], 



1 



4(1 — — y'^ — z'^) 




xy 
yz 



xz 



X 



yz 



(8) 
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obtained by the application of Dittmann's formula eq. (3.7)] (cf. 

dsuresip, P + dp)^ = ixrjdpdp + .^(dp - pdp)(dp - pAp)}. (9) 
4 \p\ 

(Of course, in the limit v ~f 1, pp, in effect, degenerates to a two-level system, and it is of interest to keep this in 
mind in examining the results presented here. In the opposite limit w — *■ 0, one simply leaves the domain of quantum 
considerations. ) 

We note that (||) differs from (^ in that it has an additional factor, 

" 4(2;2 +y-^ + z^- {v - 2)2) ' ^^^^ 

Since v can be no greater than 1 and + y'^ + z'^ no greater than w if is to meet the nonnegativity requirements 
of a density matrix, / must be negative over the feasible range (P) of parameters of pp. In fact, the square root of 
the determinant of the "complementary" 4x4 submatrix of ^ — the one associated with the nullified parameters, 
s, t, u, w, rather than v, x, y, z — is equal to /. 

Now, it is interesting to note — transforming the Cartesian coordinates (x, y, z) to spherical ones (r, 9, 4>) — that 
while the previous result ^ of |^ can be normalized to a (proper) probability distribution over P, 

^ a A\ 3r2sin6' 

p{v, r, 0, 0) = — — —r, (11) 

47r^i;(l — ti) 2 (v'^ ~ r^) 2 

the new prior (^ is itself not normalizable over P, that is, it is improper. However, we can (partially) integrate (|^) 
over the three spherical coordinates to obtain the univariate marginal over the variable v, 

/ X TT^ 2 1 , 

64t; n_y\2 -l + i; 

The integral of (p^ over v e [0, 1] diverges, however. We can compare ( p^ ) with the univariate marginal probability 
distribution of (|ll|) jl^ , eq. (19)], 

Piv) = 777^- (13) 
4(1 - v)-2 

(In | ]30| , pl| p{v) was interpreted as a density-of-states or structure function, for thermodynamic purposes, and the 
associated partition function reported.) The behaviors of ( p^ and (|l3|) are quite distinct, the latter monotonically 
increasing as v increases, while the [negative] of the former has a minimum at w « .618034. 

Let us also observe that the factor /, given in (|l^), the added presence of which leads to the non-normalizability of 
(^, takes the form in the spherical coordinates, 

^ 4{r + 2~v){r + v-2)^ 4{r^ - {v - 2)2) ■ ^^^^ 
The eight eigenvalues (A) of the nullified form of the Bures metric tensor g, that is g, come in pairs. They are 

''^1,2 = -j-i ^3,4 = — T^, A5,6 = -—, — ; rr, (15) 

4:V 4 + 2r — 2v 2[r + v — 2) 



-2(r2 + {v~ 2)v) + 2(r4 + i;* + 2r2(2 + {v ~ 4)w))5 ' 

Of course, the product of these eight eigenvalues gives us \g\, the square root of which — that is (||) — constitutes 
the new (but unnormalizable/improper) possibility here for the conditional Bures/quantum Jeffreys' prior over the 
four-dimensional convex subset P of the eight-dimensional convex set Q composed of the 3x3 density matrices. 
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II. BURES PRIORS AND SEPARABILITY PROBABILITIES FOR VARIOUS COMPOSITE QUANTUM 

SYSTEMS 



A. One-parameter 2 ® 2 systems 



Now, let us seek to extend the comparative form of analysis in sec. [B to the 4x4 density matrices. For the 
Bures metric in this setting we rely upon Proposition 1 in the recent paper of Dittmann ||^ , which presents an explicit 
formula in terms of the characteristic polynomials of the density matrices. (Let us point out that in the earlier preprint 



versions, in particular quant-ph/9911058 y4, of our paper here, a number of "anomalous" results were reported. These 
turned out to be attributable to our misinterpretation of the symbol Y' in jsj as the transpose of y, rather than the 
conjugate transpose of Y . We have since amended our analyses in this regard.) We apply it to several one-dimensional 
convex subsets of the fifteen-dimensional convex set (i?) of 4 x 4 density matrices. These subsets — unless otherwise 
indicated — are (partially) characterized by having their associated two 2x2 reduced systems described by the 
fully mixed (diagonal) density matrix, having ^ for its two diagonal entries. Or to put it equivalently, the three 
Stokes/Bloch parameters for each of the two subsystems are all zero. (A complete characterization of the inseparable 
2 (g) 2 systems with maximally disordered subsystems has been presented within the Hilbert-Schmidt space formalism 



1. The three intra- directional correlations are all equal 

For our first scenario, we stipulate zero correlation between the spins of these two reduced (fully mixed) systems in 
different directions, but identical non-zero (in general) correlation between them in the same {x, y or z) directions. We 
denote this common correlation parameter by C,. In terms of the parameterization of the coupled two-level systems 
given in Q (cf. p^ , ^ ), the feasible range of C is [^j, i^]- (The parameterization in is based on the superposition 
of sixteen 4x4 matrices — which are the pairwise direct products of the four 2x2 Pauli matrices, including among 
them, the identity matrix. Since the six Stokes/Bloch parameters have all been set to zero, the nine correlation 
parameters {Cij,hj — x,y,z) must all lie between -1 and 1, and the nine- fold sum of their squares can not exceed 
3 It has iDcen shown that all the tangent vectors corresponding to a basis of the Lie algebra — corresponding 
to two copies of SU{2) — span six dimensions, and thus there are, in fact, nine nonlocal parameters [Q. A referee 
has suggested that the use of "local orbits would simplify the picture especially for physicists dealing professionally 
with entanglement. It is because then the reader knows that one deals with what is of the main importance [orbit 
parameters — like Schmidt coefficients for pure states] from the point of view of say quantum information transmission 
[hke e. g. teleportation]" (cf. H^JH)) 

If we implement the formula of Dittmann |^, eq. (9)] using a general (fifteen-parameter) 4x4 density matrix [Q, 
then nullify twelve of the parameters of the resultant Bures metric tensor, and set the indicated remaining three {(.a) 



all equal to one value C, we obtain as the conditional Bures prior (the counterpart of \g\ 2 in sec. [B) 

32768 



(l-4C)3(l + 4C)3^/T-T2C 



(16) 



On the other hand, if we set the fifteen parameters in precisely this same fashion before employing the formula of 
Dittmann, we obtain for the volume element 



(1-8C-48C2)5 



The former prior is non-normalizable over C g l^]' while the latter is normalizable, its integral over this interval 
equalling ^. In Fig. |l|, we display this probability distribution. The pair of outcomes ((p^) and (p^)) is, thus, fully 
analogous in terms of normalizability, to what we found above ((^) and (H)) for the particular four-dimensional case 
(P) of the three-level quantum systems [Q) investigated above. 
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Bures prior 
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15 - 



10 - 



5 - 



-0.25 -0.2 -0.15 -0.1 -0.05 0.05 
FIG. 1. Normalized conditional Bures prior (^) for one-parameter four- level scenario 1 

Now, for e 1^]' associated one-parameter density matrix is separable or classically correlated (a necessary 

and sufficient condition for which for the 4x4 and 6x6 density matrices is that their partial transposes have non- 
negative eigenvalues |^^). So, if we integrate the normalized form of ( [TtI ) over this interval, we obtain the conditional 
Bures probability of separability (cf. Jl|, ^ , p] Jg|] ) . This probability turns out to be precisely i. Contrastingly, in for 
arbitrary coupled two-level systems in the fifteen-dimensional convex set i?, it was necessary to rely upon numerical 
(randomization) simulations for estimates of the Bures probability of separability, so this exact result appears quite 
novel in nature. (In pi, the [unconditional] Bures probability of separability was estimated to be ~ .1.) 



2. One intra- directional correlation equals the negative of the other two 

A closely related scenario in which the probability of separability is also precisely i is one for which the only 
non-nullified parameters are again the three intra-directional correlations, but now two of them (say, for the x and 
j/-directions) are set equal to Q and the third to —C,. Then, the conditional Bures probability distribution (computed 
in the analogous manner) is (Fig. ^ 



4V3 



7r(l + 8C- 48(2)5 



(18) 



The region of feasibility is [— t^, t] and of separability, [— t^j t^] 



Bures prior 



20 



15 



10 - 



-0.05 0.05 0.1 0.15 0.2 0.25 

FIG. 2. Normalized conditional Bures prior (nsl) for one-parameter four- level scenario 2 
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3. The six inter- directional correlations are all equal 



Let us now examine another one-parameter scenario in which the pair of two-level systems is still composed of fully 
mixed states, but for which the correlations (Qi) in the same directions are zero, while the correlations in different 
directions {Qj, i ^ j) are not necessarily zero and all equal. Thus, we ab initio set the (six) interdirectional correlations 
to C, the other nine parameters all to zero, and employ the formula of Dittmann (in the manner, we have settled 
upon for this and all subsequent analyses here). We obtain the conditional Bures probability distribution (Fig. 0), 



8V2 



7r(l - 8C - 128(2 



(19) 



over the feasible range, C G ^he range of separabihty is [— j^, t^]- The associated conditional Bures 

probability of separabihty is then ^ + — — ^ « .608173. 



Bures prior 

35 
30 
25 
20 
15 
10 
5 



-0.125-0.1-0.075-0.05-0.025 0.025 0.05 

FIG. 3. Normalized conditional Bures prior ([L9|) for one-parameter four-level scenario 3 



4- Tfie six inter- directional correlations all equal the negative of the three intra- directional correlations 

Another one-dimensional scenario of possible interest is one in which we set the intra-directional correlations to C 
and the inter-directional ones to Now, the range of feasibility is C G [~^^ J2] ^^^'^ ^^"^ interval of separability is 
^ ^ [^^' m]' ^o^j application of the Dittmann formula yields 

12 (20) 

(4C-1)(12C-1)(1 + 20C)' ^ ^ 

the square root of which gives us the unnormalizcd Bures prior. Since the integrations involved yield various elliptic 
functions, we have to resort to numerical methods to obtain the Bures probability of separability, that is, .702675. In 
Fig. 1^, we plot the associated Bures probability density function. 
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Bures prior 




FIG. 4. Normalized conditional Bures prior for one-parameter four-level scenario 4 



5. Three scenarios for which the probabilities of separability are simply 1 

If we set all nine (inter- and intra-) directional correlations to one value and the other six (Stokes/Bloch) 
parameters to zero, so that again the two reduced systems are fully mixed in nature, then proceeding along the same 
lines as above, we obtain the particularly simple conditional Bures probability distribution, 

12 . . 

(21) 



7r(l - 144(2)5 



over the feasible range C, G [— i^, til- However, all the states in this range are separable, so the associated probability 
(Bures or otherwise) of separability is simply 1. 

If we (formally, but somewhat unnaturally) set all fifteen parameters to Q, say, then the conditional Bures prior is 
proportional to 



2(3-20C)5 



(1 + 12C - 336C2 -I- 576C3 



(22) 



Though version 3 of MATHEMATICA failed (exceeding its iteration limit of 4096) to integrate over C, £ [- ^.„ \ ^. 

4(3+2 V 3 ) 

— .0386751, version 4 (as shown by M. Trott) yielded 



l\l + ^^"^A^^ + V3);sin-i _^(i3 _ 4V3))|^(13 + 4^3)). (23) 

In any case, all the 4x4 density matrices in this one-dimensional set are separable, as well. Another scenario in which 
the probability of separability is unity, is one in which the three intra-directional correlations (Cu) are all zero, and 
the two systems are aniz-correlated in different directions, that is Qj = —Qi- 



6. Rains- Smolin entangled states 

On p. 182 of |4|], Rains presents a one-parameter (x) set of 4 x 4 density matrices, apparently communicated to 
him by Smolin. The corresponding normalized Bures prior for this set of states — defined over the range of feasibility 
X e ulu = V807599/175 « 5.13523 — is 

m 



7r(807599 - 30625a;2) 
None of the members of this set is separable. 
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7. Two-qubit Werner states 



It is of some interest that all the Bures conditional probabilities of separability we obtained in the various one- 
dimensional scenarios above are substantially larger than the approximate estimate of .1 for the fifteen-dimensional set 
of 4 X 4 density matrices, obtained on the basis of (unfortunately, but perhaps unavoidably, rather crude) numerical 
methods in One does, however, obtain a (somewhat smaller) probability of separability of -j for the two-qubit 
"Werner states" These are mixtures of the fully mixed state and a maximally entangled state, with weights 1 — e 
and e, respectively. (In terms of our other set of parameters, the three intra-directional correlations are all equal to 
— e/4, while the remaining twelve parameters are zero. A referee has suggested that the parameter e is comparable to 
the visibility in optics [Q, that is {Imax — Imin) / {Imax + Imin), where / is the intensity.) The range of feasibility is 
e G [0, 1] and of separability, [0, |]. The Bures conditional probability distribution (Fig. ^) is 



3^/3 



7r(4-|-8e- 12e^)2 



(25) 



Bures prior 
6 

5 ; 

4 ^ 
3 
2 
1 



0.2 0.4 0.6 0.8 1 

FIG. 5. Bures conditional probability distribution (G5[) over the two-qubit Werner states 



B. Tuio-parameter 2 ® 2 systems 

1. Two intra-directional correlations are equal and the third one, free 



Now we modify scenario 2 of sec. II A 2, in that we set two intra-directional correlations again to a common value, 
call it and the third, not to — C this time, but to an independent parameter, call it r/. (The remaining parameters 
— the six Stokes/Bloch ones and the six inter-directional correlations stay fixed at zero.) The normalized conditional 
Bures prior is then 

(26) 



7r((l+4r?)((l- 477)2 - 64(2)) = 



The (triangular-shaped) range of feasibility over which we integrated to normalize the (conditional) Bures volume 
element extends in the 7/-direction from —jtoj. In this triangle, we integrated first over ( from ^^^'^^'^'> to -^i^^^. 
The part of the (rhombus-shaped) range of separability for rj e [0, j] coincides with the feasible domain, and for 
r] € extends over ( e [~^^^^^^, ^^^-^^^]. The univariate marginal probability distribution (Fig. ^) of ( p6| ) over 

ry is V2/^/TTW]. 
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Marginal Bures prior 



20 



15 



10 



h 

-0.2 -0.1 0.1 0.2 

FIG. 6. Univariate marginal Bures prior probability distribution for two-parameter four-level scenario 

The probability of separability for this two-parameter four-level scenario is, then, remarkably simply, \/2 — 1 « 
.414214 — being somewhat less than the ^ of the related scenario 2 of sec. [I A 2. (Of this total figure, 1-^ w .292893 

comes from the integration over r] > and ~ 2 « .12132 from the other half of the rhomboidal separability region, 
that is for 7] < 0. This second result required the use of version 4 of MATHEMATICA, and I thank Michael Trott for 
his assistance.) 



2. States inferred by the principle of maximum nonadditive (Tsallis) entropy 

Here the two variables parameterizing the 4x4 density matrices are the g-expected value {hq — "internal energy" ) 
and the variance {a^) of the Bell-Clauser-Horne-Shimony-Holt (Bell-CHSH) observable or "Hamiltonian" used 
by Abe and Rajagopal (cf. |Q) in their effort to avoid fake entanglement when only hq is employed in the Jaynes 
maximum entropy inference scheme [Q. We know from that the feasible region is determined by < 6, < 2-^2 
and 2V26g < fi^ < 8. 

Let us, first, set the positive parameter q indexing the Tsallis entropy to 1. ("It is of interest to note that for g > 1, 
indicating the subadditive feature of the Tsallis entropy, the entangled region is small and enlarges as ones goes into 
the superadditive regime, where q < 1" p^ .) Then, the corresponding Bures prior probability distribution — again 
applying the formula of Dittmann — is 



(27) 



In Fig. 1^, we show the univariate marginal probability distribution of ( p7| ) — having integrated it over the parameter 
af — for the expected value bi. The region of separability is determined | |47| , eqs. (11), (12)] by the supplementary 
requirements that af < 8 — 2\f2b\ and b\ < V2. The probability of separabihty is, then, V2 — 1 w .414214 (cf. 0, 
Fig. 1(d)]). 
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Marginal Bures distribution 
1 . 4 
1.2 



0.5 1 1.5 



2 . 5 




FIG. 7. Marginal Bures prior probability distribution over the expected value bi 

For q = i, the corresponding Bures probability distribution is (cf. (|2 

32 



Tr{32 + Abi^ + {al-8)al)2' 

2 2 2 

In Fig. ^, we show the marginal probability distribution ol (|2^) over the ^-expectation value bi 

Marginal Bures distribution 

0.8 



(28) 



0.6 



0.4- 



0.2 



b 



FIG. 



0.5 1 1.5 2 2.5 
Marginal Bures prior probability distribution over the expected value bi 



The probability of separability is then again, quite remarkably, \/2 — l. (The domain of integration is now determined 
by the supplementary requirements that cr| < 8+2\/2&i — 2\/2\/f'i (4\/2 + fei ) and 6i < 4— 2\/2.) So, it would appear 

2 2 y 2 2 2 

from our two analyses that the probability of separability is independent of the particular choice of q. (Certainly, 
whether or not this is so bears further investigation. However, we have encountered initial computational difficulties 
in obtaining results for other choices of the index q.) 

Of course, it would be of interest to consider the index of the Tsallis entropy g as a third intrinsic variable 
parameterizing the joint states of the two qubits, in addition to bq and cr^, but doing so would appear to exceed 
current computational capabilities. 



C. T/iree-parameter 2 (g) 2 systems 

1. The three intra- directional correlations are independent 

If we now modify the scenario immediately above by letting all three intra-directional correlations be independent 
of one another — while maintaining the other twelve (Stokes/Bloch and inter-directional correlation) parameters of 
the four-level systems at zero — the conditional Bures prior takes the form (symmetric in the three free parameters) 

r. (29) 

((-1 + 4C - 4?? - 4k)(1 + 4C + 4ry - 4k)(1 -I- 4C - 47? + 4k)(-1 + 4C + 4?? + 4k))^ 
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The normalization factor, by which this must be divided to yield a probability distribution is The Bures probability 
of separability is f — ^ ~ .13662. 

To obtain these results, we employed the change-of- variables, k = —j+C — rj—j, that is, v = — 1 + 4^ — 477 — 4k. The 
integration over the domain of feasibility was obtained using the ordered limits, C G [—1/4, 1/4]; v e [2(— 1 + 4(^),0]; 
and T] £ [(-2 - v)/8, {&( - v)/8]. For the domain of separability, we employed C e [0, 1/4]; v e [2(-l + 4C), 0]; and 
77 e [(—2 + 8( — v)/8, ~v/8]. The univariate marginal probabihty distribution of ( |29| ) over ( G [—1/4, 1/4] is the 
uniform one. 



2. Diagonal density matrices 

As a simple exercise of interest, we have analyzed — again with the use of Dittmann's general formula for the 
Bures metric tensor — four-level diagonal density matrices, having entries denoted x,y, z,l — x — y — z. (All such 
density matrices are separable.) The Bures prior probability distribution over the three-dimensional simplex spanned 
by these entries is, then, simply the Dirichlet distribution 



Tr'^(xyz{l ^ X — y — z)) 2 

which serves as the prior probability distribution ("Jeffreys' prior") based on the Fisher information metric for a 
quadriinomial distribution p^jl9| ]. This result is, thus, consistent with the study of Braunstein and Caves [^, in 
which the Bures metric was obtained by maximizing the Fisher information over all quantum measurements, not just 
ones described by one-dimensional orthogonal projectors. 



D. A Fojtr-parameter 2 (g) 2 system 

1. One-parameter unitary transformation of diagonal density matrices 

If we transform the (three-parameter) diagonal density matrices discussed immediately above by a (one-parameter) 
4x4 unitary matrix, 

1 0^ 



„P I COS w sin w 
— sinui cosw 
.0 1, 



C/i=e™^= " "''"^ " (31) 

' — smui cosw ' 



where P is a traceless Hermitian matrix (one of the standard generators of S'[/(4) [^) with its only nonzero entries 
being i in the (3,2) cell and —i in the (2,3) cell, we find that the (unnormalizcd) Bures prior is (cf. (|30|)) 



Viv^ (32) 



8{xyz{y + z)(l - x-y- z)) : 



being independent of the fourth (unitary) parameter w. To normalize this prior over the product of the three- 
dimensional simplex and the interval w € [0, 2tt], we need to multiply it by The additional separability requirement 
is that w G [— w,m], where 



sm (-^^7 = ) 

u= ^ . (33) 

The approximate Bures probability of separability is, then, .112, similar to the (unrestricted) estimate in [^. 
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E. One-parameter 3 (8i 3 systems 



— the two- qutrit Werner states 



Caves and Milburn |5^, sec. Ill] have constructed two- qutrit Werner states. (Such states violate the partial 
transposition criterion for separability, while satisfying a certain reduction criterion, the violation of which implies 
distillability ^^■) Again applying Proposition 1 of Dittman ||] to this one-dimensional set of 9 x 9 density matrices, 
we obtained for the (unnormalized) conditional Bures prior, the square root of the ratio of 

- 16(2 + 7e) (496 + 14384e + 179472e2 + 1269568e^ + 5676488e'^+ (34) 



to 



16753596e^ + 314196466*^ + 31863023e^ + 14859999e 



3(-l + ef{l + 8e)(31 + 161e)(31 + 603e + 39936^ + 89816^ 



The range of separability is [0, i] [Q. (If the maximally entangled component of the Werner state is replaced by 
an arbitrary 9x9 density matrix, the range of separability for the resulting mixture must include [0, ^] ||5C|].) The 
integral of the square root of the ratio (||) over this range is 1.05879, while over [0, .999999], it is 9.62137 x 10^. The 
conditional Bures probability of separability of the two-qutrit Werner states, thus, appears to be vanishingly small 
(cf. 1^,0. The "maximally entangled states of two qutrits are more entangled than maximally entangled states of 
two qubits" [5C| ]. The two-qutrit Bures prior (Fig. ^) is much more steeply rising than the two-qubit one displayed in 
Fig. |. (Note, of course, the difference in the scale employed in the two plots.) 



Bures prior 
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0.999 0.9992 0.9994 0.9996 0.9998 
FIG. 9. Bures conditional (unnormalized) measure — that is, the square root of the ratio 
states 



over the two-qutrit Werner 



F. One-parameter 2 ® 3 systems 

It clearly constitutes a challenging task to extend (cf. sec. p^I B| ) the series of one-dimensional analyses in sec. || 
to m-dimensional (m > 1) subsets of the fifteen-dimensional set of 4 x 4 density matrices, and a fortiori to n x n 
density matrices, n > 4. Only in highly special cases, does it appear that the use of exact integration methods, such 
as exploited above, will succeed, and recourse will have to be had to numerical techniques, such as were advanced in 
. We should note, though, that in those two studies, numerical integration procedures were not readily applicable. 
They would be available if one has, as here, explicit forms for the Bures prior, and can suitably define the limits of 
integration — that is, the boundaries of the sets of density matrices (and separable density matrices) under analysis. 



1. Scenario 1 



As an illustration of the application of numerical integration techniques to such a higher-dimensional scenario, we 
have considered the one-parameter (y) family of 6 x 6 density matrices. 
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(35) 



the 2x2 and 3x3 reduced systems of which are fully mixed. The range of feasibility \s v <^ [—.0546647, .10277] 
and that of separability is G [—.0546647, .0546647]. (Again, we apply the positive partial transposition Peres- 
Horodecki condition, sufficient as well as necessary for both the 2 (g) 2 and 2 (g) 3 systems [^.) We have determined 
the corresponding (conditional) Bures prior, again based on Proposition 1 in [^. The probability of separability is 
the ratio of the integrals of the prior over these two intervals. This turns out to be .607921 — as seems plausible from 
the plot of the (conditional) Bures prior over the feasibility range in Fig. |l^. This figure displays a minimum at = 
(cf. Fig. For < 0, we numerically integrate, as well as plot, the absolute value of the (imaginary) Bures prior in 
this part of the parameter range. 



3ures prior 
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0.025 0.05 0.075 



. 1 



FIG. 10. Normalized conditional Bures prior over the single parameter (y) of the sis-level system (|35| 



2. The qubit-qutrit Werner states 



The probability of separability also appears to be vanishingly small for the "hybrid" qubit-qutrit Werner states (cf. 
[^). The greatest degree of entanglement one can hope to achieve is to have the qubit in a fully mixed state, and the 
qutrit in a degenerate state, with spectrum ^, ^ and 0. The range of separability can then be shown to be 0, j] [^ . 
(For the analysis of multi-quhit Werner states, see [^.) The (unnormalized) conditional Bures prior (Fig. y) is the 
square root of the ratio of 



10(1 + 2e)(26 + 286e + 1236e^ + 2506e^ + 2021e^) 



(36) 



to 



(-1 + e)2(l + 5e)(13 + 32e)(13 + 126e + 4296^ + 5l2e^). 
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III. CONCLUDING REMARKS 



We believe the main contribution of this study is that it emphatically reveals — through its remarkably simple 
results of an exact nature — the existence of an intimate connection between two major (but, heretofore, rather 
distinct) areas of study in quantum physics. By this is meant the study of: (1) metric structures on quantum 
systems Q; and (2) entanglement of quantum systems (cf. @). It should be noted, however, that the several 
exact Bures probabilities of separability adduced above are all for certain qubit-qubit systems (representable by 4 x 4 
density matrices). It would be of interest to see if exact results are obtainable for qubit-qutrit (representable by 
6x6 density matrices) and even larger-sized systems, as well as for qubit-qubit systems parameterized by more than 
three variables (the maximum possible being fifteen). Such investigations would, undoubtedly, demand considerable 
computational resources and sophistication. (Let us also direct the reader's attention to our study |9|, in which we 
employ numerical methods to estimate the Bures probability of separability of the two-party Gaussian states, forms of 
continuous variable systems. For this purpose, we employed an analogue |^ (cf. |Q) of the Peres- Horodecki criterion 
for separabihty [^ .) 

We also intend to use the Bures probabilities as weights for the entanglement of formation [^6| . By doing so, we 
should be able to order different scenarios by the total amount of entanglement they involve. As a first example, we 
have found this figure to be equal to .0441763 for the first (one-parameter) scenario we have analyzed above (sec. [I A), 
that of three equal intra-directional correlations (and all other parameters fixed at zero). 

One of the clearer findings of 0| was that the Bures (minimal monotone) probability of separability provides an 
upper bound on the related probability for any monotone metric. So, it would appear that the results reported above 
might also be interpreted as providing upper bounds on any acceptable measure of the probability of separability. 

Let us also remark that a quite distinct direction perhaps worthy of exploration is the use of Bures priors as 
densities-of-states or structure functions for thermodynamic purposes (cf. pO|j3l| , |59| -|6T[ ) . (However, our initial efforts 
to find explicit forms for the partition functions corresponding to the results above have not succeeded.) 
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